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TECHNICAL NOTE 3134 


A METHOD FOR ESTIMATING VARIATIONS IN THE ROOTS OF THE 
LATERAL-STABILITY QUARTIC DUE TO CHANGES IN MASS 
AND AERODYNAMIC PARAMETERS OF AN AIRPLANE 
By Ordway B. Gates, Jr., and C. H. Woodllng ' 

SUMMARY 


A method is presented for estimating variations in the roots of 
the lateral-stability quartic due to changes in mass and aerodynamic 
parameters of an airplane. The method is applied to three high-speed 
airplanes and the changes of their lateral stability characteristics 
are determined by considering increments in various airplane parameters. 

The expressions indicating the rate of change of the Dutch roll 
damping and frequency with respect to the stability derivatives and 
mass characteristics are simplified and the results when compared with 
those obtained from the exact expressions show very good agreement. 

The rates of change of the roots with reBpect to the parameters 
are shown to have a definite relationship with the amplitude coeffi- 
cients and ratios of the lateral modes of motion of the airplane sub- 
sequent to applied forces or moment s . From these relationships, calcu- 
lation of the rate of change of the roots with respect to five prescribed 
parameters allows determination of the remaining partial derivatives, 
the amplitude coefficients, and ratios mentioned above. 

The derived expressions should afford some insight into the types 
of automatic stabilization devices likely to be most effective for a 
given airplane since. If automatic-stabilizer dynamics are neglected, 
the stabilizer is effectively varying one or more of the mass or aero- 
dynamic parameters of the airplane. 

A method is given in the appendix which can be used to calculate 
approximately the roots of the lateral- stability quartic. 
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UJTROnJCTION 


Recent investigations of the lateral stability characteristics of 
airplanes have indicated that small variations in the esti m ate of the 
mass and aerodynamic parameters of a given airplane may cause pronounced 
changes in its stability. (For example, see refs. 1 to 4.) Estimations 
of the mass and aerodynamic parameters of the airplane, "whether from 
wind-tunnel data, flight tests, or existing theory, are subject to 
certain probable errors. Hence, a means of evaluating the effect of 
these anticipated probable errors on the stability characteristics of 
a given airplane should prove very useful. Also, such a tool should 
provide some insight into the relative importance of parameters or com- 
binations of parameters affecting the stability of the airplane and, as 
a result, provides trends which should be useful in the selection of 
automatic stabilizing devices for particular airplanes. 

The purpose of the present investigation is to derive expressions 
from which the approximate variation of the roots of the lateral- 
stability characteristic-quartic equation due to small changes in one 
or more of the airplane parameters can be calculated. As a check, the 
results obtained by utilizing the derived expressions are compared with 
the results of exact calculations of the variation of the roots due to 
changes in certain parameters for three different airplanes. Also, 
there is presented a method which can be used to calculate approximately 
the roots of the characteristic-quart ic equation for a given set of 
parameters. 

The rates of change of the roots of the characteristic stability 
quartic with respect to the airplane mass and aerodynamic parameters 
are shown to have a unique relationship with the amplitude coefficients 
and the ratios of the modes of the lateral motions of the airplane sub- 
sequent to certain disturbances. 

K. Mitchell in reference 5 presented a method for calculating the 
approximate changes in the roots of a quartic equation due to variations 
in parameters upon which its coefficients depend which is very si m ilar 
to the present method. The present investigation, conducted independ- 
ently of that of Mitchell, makes use of a slightly different approach 
to the problem and, in addition, the derived expressions are examined 
with a view toward obtaining simpler expressions with which to work. 
Also, an attempt is made to establish a criterion which defines the 
range of variations in the parameters for which the assumption of 
linearity is valid. 
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SYMBOIS AMD COEFFICIENTS 


t 


A,B,C,D,E coefficients of lateral- stability equation 

B ' ,C' ,D' ,E 1 coefficients of lateral- stability equation having been 
divided by A 

a measure of damping of Dutch roll oscillation 

a ± iru complex root of stability equation corresponding to the 

Dutch roll oscillation 


a,b,c,...i coefficients in solutions of lateral equations of motion 


b 


■wing span, ft 


c L 


trim lift coefficient, 


W cos y 


Cl 


rolling-moment coefficient. 


Rolling moment 
qSb 



°*r 



Cn 


dancing- in- roll derivative, rate of change of rolling- 

§ oment coefficient with rolling-angular-velocity factor, 
C, 

— per radian 

^gb 

2V 


rate of change of rolling-moment coefficient with yawing- 

307 

angular-velocity factor, — pg-, per radian 

S— 

2V 


effective-dihedral derivative, rate of change of rolling- 
moment coefficient with angle of sideslip, dCj/Sp, per 
radian 


yawing-moment coefficient. 


Yawing moment 
qSb 



rate of change of yawing-moment coefficient with rolling- 


3Cn 


per radian 


2V 


angular- velocity factor, 
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Cn r 


Cnp 


C Y 

Cy 

L P 

cy r 


% 

Db 

F(A) 

g 

i 


K Xo 


K xz 


K Z 


% 


damping-in-yaw derivative, rate of change of yawing-moment 

BCv, 

coefficient with yawing-angular- velocity factor, — 
per radian m 


directional-stability derivative, rate of change of yawing- 
moment coefficient with angle of sideslip, dCn/dp* per 
radian 


lateral- force coefficient. 


Lateral force 
qS 


rate of change of lateral- force coefficient with rolling- 

3Cy 

angular-velocity factor, — *■, per radian 

°2V 


rate of change of lateral-force coefficient with yawing- 
angular- velocity factor, — per radian 

% 


lateral- force derivative, rate of change of lateral-force 
coefficient with angle of sideslip, dCy/dp, per radian 

differential operator, d/dB^ 

characteristic-quartic equation 

acceleration due to gravity, ft/sec/sec 


nondlmensional radiu s of gyration in roll about longitudinal 
stability axis, ^Kx o 2 cos 2 t] + K>?; o 2 sin 2 T] 

nondlmensional radius of gyration in roll about principal 
longitudinal axis, fc^/b 

nondlmensional product-of-inertia parameter, 

(Kz o 2 " K Xo 2 )sin T] cos tj 

nondlmensional radi us of gyration in yaw about vertical 
stability axis, ^cos 2 ^ + Kx Q ^sin 2 Ti 

nondlmensional radius of gyration in yaw about principal 
vertical axis, k£ 0 /b 

radius of gyration in roll about principal longitudinal 
axis , ft 
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m 


P 

Pl>P2 


P 

q. 

< 3 - 1 j<1 2 


r 

S 

8 b 

Tl/2 

t 

V 

v 

W 

*1 

P 

7 

n 


radius of gyration in yaw about principal vertical axis, ft 


mass of airplane, W/g, slugs 

rolling velocity, &0/dt, r adlans/sec 

coefficients of A in two quadratics that constitute 
stability quartic 


Laplace transform variable, f(p) = / e“ pt F(t)dt 

J 0 


dynamic pressure. 


PV 2 

~TT’ 


lb/sq ft 


constant coefficients in two quadratics that constitute 
stability quartic 

yawing velocity, dy/dt, radians/sec 
wing area, sq ft 

nondimens ional time parameter based on span, Vt/b 

time for amplitude of oscillation to change by a factor of 2 
(positive value indicates a decrease to half- amplitude, 
negative value indicates an increase to double amplitude) 

time, sec 

airspeed, ft/sec 

sideslip velocity along lateral axis, ft/sec 
weight of airplane, lb 

mass and aerodynamic parameters upon which A,B,C,D, and E 
depend 

angle of sideslip, v/V, radians 

angle of flight path to horizontal axis, positive in a 
climb, deg 

inclination of principal longitudinal axis of airplane with 
respect to flight path, positive when principal axis is 
above flight path at nose, deg 



6 


NACA TN 5134 


A,A n nondimens ional root of stability equation, 

AA^ + BA^ + CA 2 + DA + E = 0 


Ai real root of stability equation corresponding to spiral mode 

Ag real root of stability equation corresponding to damping-in- 

roll mode 

Pb relative-density factor, m/pSb 

p mass density of air, slugs/cu ft 

0 angle of bank, radians 

f angle of yaw, radians 

cjo frequency of Dutch roll oscillation 


cqn 


natural frequency of Dutch roll oscillation, (a 2 + cu 2 ) 1 / 2 

frequency of Dutch roll oscillation based on time in 
seconds, <x>' = radians/sec 


The subscript o is used to indicate the value of the quantity 
when Ax.£ = 0 . 


ANALYSIS 

Equations of Motion 


The lateral nondimensional linearized airplane equations of motion 
for level flight (7 = 0 ) and controls fixed referred to the stability 
axes are 

Rolling: 


+ Kxz 3 ^ 2 ^) = + \ C 2 p I> b^ + | 

Yawing: 

2pb(% 2 Db 2 * + = CnpP + \ 0^0 + \ C^t 
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Sideslipping: 


^(^bP + = °Y p P + J Cy p Db0 + § Cy r I>bt + Cl0 

When 0 o e^^ is substituted for 0, TjfQe^ 3 for if, and [Joe^ 3 

for 0 in the equations written in determinant form, X must be a 
root of the stability equation, 

F(A) = AK k + BAS + CA 2 + DX + E = 0 (l) 

where 

A =.8^ 3 (k x 2 Kz 2 - K^ 2 ) 

B = “ 2 ^b 2 ( 2K X 2K Z 2c Yp + Kx^nj. + - SKxz 2 ^ ~ Kxz c Z r “ KxzCn^ 

C = H b (K x 2c nr C Yp + i^-bK x 2 C np + Kz^jpCyp + | C nr C lp - - 

^b K XZ C Zp - c rLp K XZ C Y p + K XZ C np C Y p " Kx^Y^np + K XZ c Y r C Z p “ 

ls%-*p%) L 

D “ - r C %°lp% - «>%% + f c np c l r % + + 

^AzS - ^ l% 2Cj p + i o Ip C np o Yr - i C^Oj^ - 
i C !rS C Ip * i 

S ‘ I ^(V^p - °* r %) 
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Derivation of Expressions for the Rate of Change of the Roots 


With Changes in Airplane Parameters 


The lateral-stability quartic is given by equation (l). The coef- 
ficients A, B, C, D, and E have been shown to be functions of the 
various mass and aerodynamic parameters of an airplane, henceforth 
designated as xq, X£, X3, . . . Xq. If the parameters xq are assumed 

to be independent variables, a change in a parameter will result in 
changes in the coefficients of equation (l), that is. 


A - Aq + — Axq 

OXj_ 


_ „ 8b a 

B = B 0 t - — Axq 

OXq 


c = c 0 + Axi y 

oxq 




E = E 0 + — Axq 
0 8X ± 


(3) 


The roots of equation (l) will then become A = 'Kq + AA. For Axq = 0, 
the increment in the root AA is equal to zero and A = Aq, B = B 0 , 

C = Co, D = D 0 , and E =» E 0 . Substitution of equations (3) for A, B, 

C, D, and E and A * Aq + AA into equation (l) and consideration of 
the limiting case for which Axq and AA approach zero can be shown 
to be equivalent to differentiating equation (l) implicitly with respect 
to xq. Therefore, 



^■AqAo^ + 3BoAq^ + 2CqAq + D 0 


GO 
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A point of interest is that equations 
in the coefficients regardless of the 
are linear functions of the parameter 
parameters except Kxy,. 

Equation (4) expresses the slope or rate of change of the roots of 
equation (l) with respect to any parameter x^. If several parameters 
are varied simultaneously, the total change in A is given hy the 
expression 


(3) represent the exact changes 
size of Axi if the coefficients 
Xj_, which is true for all the 


dA = 






+ dx.. 




(5) 


or 


AA 


dA 

dx n 


Axj_ + 


dA 

Sxc 


Ax2 




General Foim of Variation Expressions and Criterion for 
Determining Approximate Range of Linearity 

If the coefficients of the stability quart ic, the roots of this 
quart Ic, and the partials of the coefficients with respect to a selected 
parameter are known, the rate of change of any root with respect to this 
parameter can he found from equation (h) . The values of the coefficients 
for a particular airplane can be found hy evaluating equations (2) . The 
roots may be found by either solving the characteristic quartic by var- 
ious conventional methods (for example, ref. 6) or by using the approxi- 
mation method presented in the appendix. The general forms of the 
partial of the coefficients with respect to several of the parameters 
x-l are given In table I. The parameters that are used are Cjp, Cj r , 

c Z p > c n p > c n r > c n^ C Y^> *1 > K X 0 2 > %Z 0 2 ' parameters K^ 2 and 

K^ q 2 (the principal radii of gyration squared) were chosen rather than 

Kx^ and Kg 2 since errors in Ky 2 and Kg 2 could result from errors 
in either K Xq 2 , K^ 2 , or tj. 

After the parameter Is chosen and the values of the coefficients 
and the values of the partials of these coefficients are substituted, 
equation (4) becomes a function of A. Now depending on the nature of 
the A of interest, equation (4) may be real or complex. 
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For X real, the approximate change of X due to a change in one 
of the above parameters is 


AX » Axi (6) 

ax^ 

and the new root due to the Axi is approximately 


X « Xo + |^- Axi ( 7 ) 

OXj_ 

For X complex, that is, a ± im, equation (k) becomes congplex, and 
by equating real to real and imaginary to imaginary, the condition 
exists that 


dX _ 8a_ + ± ch> 
Sxi bx-i Sxi 


( 8 ) 


and the approximate change of a and cu due to a change in one of the 
parameters is 


Aa «* ^- a Ax. 
ox-? x 


Ad » Axi 
dx^ 


Then, the new complex root due to the Ax-j_ is approximately 


x ac + g- AX! + i(^ + ^- Axi) 


(9a) 

(9b) 


( 10 ) 


It may be well to mention here that equation (It-) evaluated for 
X *= a - ico will merely give the conjugate result of equation (it-) 
evaluated for X = a + icu, and the change in -o> may be expressed as 
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or 


Ad ~ Axi 

OX^ 


which is identical to equation ( 9 b). Therefore, in the remaining dis- 
cussion when the complex root or its slope Is discussed, the complex 
root Is assumed to he a + ion 


Naturally, the slope expression will not accurately approximate 
the change of A for large values of Ax-j_. Therefore, a criterion 
must he established that will indicate the value of Axi at which the 


expression 


+ 


dA 

SET 


Axj_ 


no longer gives a suitable approximation to 


the new root. 


Since A is a function of X£, A(xi + Axi) may he obtained from 
the Taylor's series expansion in the vicinity of x^. The first three 
terms of the expansion are 

A(xi + Axi) = A(xi) + V(xi)Axi + ^ (xi)(Axi)_ + # . , (ll) 

2 * 


where A 1 , A", and so forth, denote differentiation of A with respect 
to xp In the previous part of this paper the third term and all 
succeeding terms in the above expansion were assumed negligible. Now 
it is assumed that, whenever the third term becomes sufficiently large 
to be 10 percent of the sum of the first two terms, the initial assump- 
tion of linearity no longer presents a suitable approximation to the 
new root . This criterion may be expressed as 




+ 

OXj_ 


< 0.10 


( 12 ) 


The second partial of A with respect to xi may be obtained by 
differentiation of equation (4). 
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d^A_ 

dx ± 2 


-2 ^-j(6AA 2 + 3BA + C)^- + 4a x A 5 + 3a 2 X 2 + 
4AA5 + 3BA 2 + 2CA 4 - D 



When A = a 4- io>, the above expression becomes a complex number and 

and — are equal to the real and imaginary parts of this number, 
dx ^ 2 



respectively. In order to apply the criterion (eq. (12)) for A = a 4 - icu, 
the terms a and cd must be considered separately; that is, for a change 
in cd equation (12) becomes 


1 d3n /. \2 

2 d^^ 

®b + 


< 0.10 


Since equation (l2) employs only three terms of equation (ll) it is 
possible that the Ax^ calculated from this expression will actually 
give slightly more or less than 10 percent. Although no results are 
presented, this criterion was found to be adequate for the prediction 
of nonlinearity of the curves. 


Simplification of the Slope Equation for the Oscillatory Root 

Although variations in the mass and aerodynamic parameters of an 
airplane may cause changes in all the roots of the stability quartic, 
usually changes in the oscillatory root are of pr im a r y interest; that 
is, the change in the damping and frequency of the Dutch roll oscilla- 
tion of the airplane is of most concern. Therefore, an attempt is made 
to simplify equation (4), for the case where Aq is equal to the 
oscillatory root (a * icu) , in order to reduce the amount of work required 
to obtain values for da /&x± and dto/dx^ and also to provide some 
insight into the parameters which affect these partial derivatives. 

The denominator of equation (4) is merely the derivative of equa- 
tion (l) with respect to A. If the roots of equation (l) are, as in 
the usual caBe, two real roots, namely, A^_ and Tkq, and a pair of 
conjugate complex roots, a ± ia>, the denominator of equation (4) may 
be written as 
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p'(X) = Aq j(A 2 - 2 aX + a 2 + o) 2 )(A - X-^ + (X 2 - 2 aX + 

a 2 + a) 2 ) (X - Xg) + (X - Xi)(X - X2X2X - 2a^j ( 13 ) 

Nov, If X s a + Itu is substituted into equation (l 3 )> 

F'(a- + 3 xd) = Aq |Oel + io) - X]_)(a + i^> - X 2 )( 2 o)if| 

as vas pointed out in reference 5 • 

Also, if X = a + io), the general form of equation ( 4 ) is 

5 (a + io) _ d ■+ ci 
&Xj_ e + fi 

By rationalizing the above equation and by using equation ( 8 ), 

8a _ de + cf (15) 

8 xi e 2 + f 2 

8o ce - df (16) 

8x ± e 2 + f 2 

Ohe exact form of d and c from equation (h) are 

d . . ^ - 6AP ♦ » - - 3- 2 > - l^ 2 - ^ 


C = - |^-(^a 5o) - 4ao)5) - M_(3a 2 o) - cd?) - t^-(2sm) - 
oxi dx. ± ayi ± 
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The form of e and f from equation (l 4 ) are 


e 



f 


= 2Ao0j -a(Ai + A2) 


+ a' 


2 - 0)2 



The simplifications are now made that a and are approximately 

zero compared with A2 and 0; that is, the damping of the Dutch roll 
oscillation and spiral mode are negligible compared with the damping in 
roll and frequency. Hence, d, c, e, and f become 

a = - 

dx* dx* Sx± 


c 


5 b 

5 x^ 



0 ) 


e = 2 AoO)?A 2 
f = - 2 Ac v? 


Then, equations (15) and (16) may be written as 


5 a 

dxi 



L 5c 2 

+ (XT - 

5xj_ 




2Aoa)^(A2^ + a)2) 



(IT) 
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&L-) + (- ^ 0 > + |S_ 0J 2 

ox J \ dXi ox-l 

2Aoii(X2^ + a* 2 ) 


Unfortunately, because of differences in sign of first-order terms 
•which make second-order terms important, expressions (17) and (l8) will 
not hold for a few of the parameters . Separate expressions are presented 
for these parameters in which some second-order terms (involving a and 
powers of a) are retained. 

It was found that, for the three airplanes used in this paper, equa- 
tion (17) may be used to calculate the rate of change of a with reBpect 
to Cn r , Cj r , Cnp, C^, Cy p , and t } with good results. 


For Cjp use 


V^P 30 hJj [fop 


3a 2 ) - 


dD ac 


ac, ac, 

L V 


2A c j3d 2 (X2 2 + to 2 ) 


2 2 

For and use 


2a_ = 
axi 


[gysr 1 


(6a 2 - cd 2 ) + fo-(3a) - 2af^— (6a 2 + oc£) + 3a 

a^ J [a Xl 


as ap _ Sib ( 

Sixi axi 3x^ 


2Aq(A 2 2 + CD 2 ) 


Equation (l8) may be used to calculate the rate of change of o> with 
respect to C^, C np , Cj r , C j , and ij. 

For C n and 0y R use 


X 2 ^-(co2 _ 3a2) - ^-(2a) - g- 

[8xj_ raq. 

2A^(7\ 2 2 + n 2 ) 




a? - 5-' 

dx i/ 



1 6 


MCA TN 3134 


For C 2 and Kx ^ use 


is_ . 
0 Xi 


X 2 


^-(W) + #-{ 0,2 . 3a 2) - |C_(2a) 


fol 




&X£ 


- Irl + + |H3“> 2 ) + ir-<“ s ) - 

aacjJ [dxj _.- QXj dx* dXj J 


2*033(5^^ + 


(22) 


The preceding expressions were simplified under the assumptions that 
Ap and a are small and hence are not expected to give a good approxi- 
mation to the exact values for airplanes where the spiral root Aq is 
of the same order as the damping-in-roll root A2 or for airplanes 
which possess a large degree of damping in the Dutch roll oscillation 
(large -a/con) . For an airplane where these conditions exist the 
general slope equation (eq. (4)) should he used. 


Other Possible Simplifications for Cn r and C n p 

Many times , by closer examination of the separate terms of equa- 
tion (17) set up for a specific parameter, an extremely simple equation 
may be written down which very well approximates the results obtained 
by the complete slope equation (eq. (4 ) ) . Two such examples using C nr 
and Cn p are presented. 

If the partial derivatives of the coefficients with respect to C nr 

(table I) are substituted into equation (l/), the following equation Is 
obtained: 


Sc 


+ £ %> 2 - i cl ° J ■ 0)2 ((-^b^vy + y ip % - c ifJ c Yp ) 




2A<ju^(X 2^ + ui?) 


( 2 ?) 


Two assumptions are made: First, the term Kx^Cyq is negligible com- 

pared with - and, second, the term C^Cyp is negligible compared 

with Cj p Cy p . 
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If the following two conditions are satisfiedj 



and 


Is! > 


2 • 5^*7, Ksrp 


and if Aq and X2 can he well-approximated hy 


Aq ~ 8p. b 3K x 2 K2 2 


*2 





^h% 2 


equation ( 23 ) reduces to 


da 1 

■ ■ « 

8p b K2 2 


(24) 


This result is the same as would he obtained hy assuming one degree 
of freedom in yaw which may he expressed as 


(2p b K z 2 D 2 - | C^D + = 0 


where the real part of the complex root is 



8ft)Kz 2 


a 
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and if this equation is differentiated with respect to C nr 

8a _ 1 

^ n r Spb 1 ^ 2 


Equation (17) set up for Cnp hy using the partials of the coef- 
ficients from table I gives 



x 2 

j^b( K XZ C Y p + \ C l r ) 

Ito 2 

- CO 2 j 

4u b 2 %2p 2 - 

' £ C l r C Y p - rA p + £ C l? c Y r ) 

2 Ao^(x 2^ + to 2 ) 

1 


(25) 


Again, an assumption is made as to the order of magnitude of several 
terms. The terms i CjpCy r and - i C^Cyp are assumed to be negligible 
compared with 

If the following condition is satisfied: 

2,5C lp Kz 2 (KxzCy p + | C z ^ + 10p b K x 2 KxzC np 

^(K X 2 K2 2 + Ky^ 2 ^ 



expression (25) may be well-approximated by 


8a „ 


(26) 


As has been pointed out before, these two expressions were derived 
under the assumption that certain conditions must be satisfied and, if 
these conditions are not satisfied for a particular configuration, equa- 
tions (23) and (25) should be used. 

It is interesting to note that, from equation ( , 2 k ), the effect of 
the derivative Cn r is highly dependent upon the moment of inertia in 
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yaw. Also, from equation ( 2 . 6 ) the derivative C n appears to be most 

P 

effective for airplanes with a low moment of inertia in roll or a high 
effective dihedral Cj ^ . 


Relation of Slope Equation to Amplitude Coefficients 
of Lateral Modes of Motion 


The Laplace transformation (ref. 7 ) of the lateral equations of 
motion for zero initial conditions and the assumption of an input C n = 1 
yields the following equations: 


^2pb%*T? " ^ C nr ^pi^(p) + ^ Cj^p^Cp) - CnppCp) - ~ 


(Wb K XZP - i C* r )pt(p) + (V b K x 2 p 2 - | C lp p)0(p) - Cj p 3(p) = 0 
(2Pb ~ | c Y r )pt(p) + (-CL - | C Yp p)0(p) + (2nbP - Cy p ^)p(p) = 0 


> 




(27) 


where p is the Laplace transform variable. 

By use of determinants, expressions are obtained for pijr(p), 0(p) , 
and p(p) as follows: 


Pt(p) 


Pl(p) 

pQ(p) 


0 (p) 


p 2 (p) y 
pq(p) I 


p(p) 


p 3 (p) 

pQ(p) 


(28) 
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where P-|_, P 2 , Pj, and Q are functions of the mass and aerodynamic 

characteristics of the airplane. By use of the Heaviside expansion 
theorem (ref. 7) the inverse transformation of equations (28) yields 
the solutions Dijr, f, and j3 as functions of the nondimensional time 
parameter s-^, that is, 


Pi(o) P^) Vb 

DtJt = — + 2_ — =— = — e - 

Q(o) n=L 7s.n0.* (An) 

A^S^ AjflS-);) 

D>jf = ao + a^e + . . . +■ a^e 


0 = .. V“ p 2(A n ) J'n sh 


Q(o) n-1 ^'(An) 


f - t> 0 + 


AiSb 


+ . . . + hme 


Am s b 


( 29 ) 


P 3^°) + v 3 ~ ^(Ap.) e A n Sb 

Q(°) n=l A n Q (A n ) 


3 


+ c 


+ 


+ 

■> 


Hie linear and distinct roots of the lateral-stability quartic 
Q(p) = 0 are A^, Ag, . . . A^ and Q, 1 denotes differentiation with 

respect to p. The roots of Q(p) = 0 are, of course, identical with 
the roots discussed previously for equation (l). 
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The following relationships can be shown to exist : 



(30) 


From the expression for the conclusion can be made that, for 

airplanes with a high 0/ \Jc ratio in the Dutch roll, the Dutch roll 
damping will be more sensitive to the derivative C^p than to the 

derivative C nr . 

Comparable relationships can be derived for the derivatives Ci r , 
C V Cjp, and Cy^ for inputs Cj = 1 and Cy = 1. For C j = 1, 
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and therefore 


m 


D\jr = do + 2 _ d n e 
n=l 


^n s b 


0 = e D + ^ e n e XnSb f 


n=l 


X s - AnS-fc 

P “ f o + 2^ f n e 


_ 1 , , 

5^-5 W. 


d\i _ ^n 2e n 

* 1 , 


P 




— ‘V-^n 


> 


r 0\ „ 2. 

\ P /A=A n 'Vl dAn/dCjg 

W = ^/^l p 

^/A*A n 

t> - \ 

2 a\ 1 /ao lr 


(31) 


(32) 
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For Cy = 1 , 


and therefore 



n=l 


0 = ho 


m 


h n e 


^n B b 


> 


rhjr = 



^ 





a a W^ P 

SA n /dCYp 




/0\ s ^/dC Y p 
WA=A n " ^/SCY r 

/g \ s Al 5A n/^Y B 
2 dAa/dtJy^ 
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Kie coefficients Eq, h 0 , c d , d 0 , e Q , f o; g Q , h D , and i 0 which 
appear in equations ( 29 ), (3^), and- (33) are defined as follows: 
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The following relationships exist for 


Kx 2 , And Kxz;: 


<&n 

SKz 2 


“Vb^n 


^n r 


dK X 2 


^b\i 


5 \i 


> 


SKxz 


-^M-b^nl 




(3 6) 


These partial derivatives were not discussed previously hut are related 
to — — , and as follows: 

^o 2 S% 0 2 


dA n _ d K Z 2 + 5}\ a dK x 2 + ^%2; 

SK Zq 2 SK Z 2 SK2 q 2 dK X 2 dK^ 2 d%2 Sk Z q 2 


= ^Kg g + ^ Sk X 2 + ^ K XZ 

^k Xq 2 Sk Xo 2 Sk x 2 Sk Xq 2 aK xz Sk Xo 2 


> 


(37) 


^ = fon_ ^n_ d%Z 

^ Sk z 2 Sk x 2 ^ ^ K xz ^ 



2 6 


NACA TN 313^ 


Substitution of tbe following relationships: 


-K - 

s COS 2 T] 

SK X 2 
— ^-15 

= sin 2 r] 

7 ? 

= sin T] cos tj 



dK z 2 
n 0 


a %0 


OKz 2 

7\ ~ 

s sin 2 r) 


= COS 2 T) 

8Kxz 

= -sin r] cos rj 

& 

iff 

0 

r\ 


aK X 0 


a K Xo 2 


dK z 2 _ 

“2Kxz 

3k x 2 _ 
dtl 

2KXZ 

^ K xz 

•I 

iff 

K> 

1 

iff 

to 


and those of equations (3 6 ) into equations (37) yields 


dA- 




SAn 





dA 




cos 2 


" + sln2 ' 1 + (t^ + ir) 8in " 008 11 


dK 


-^b A n 


*o 


— ^ — sin^T) + cos 2 tj - -\sln q cos ql 

K W *v _1 


- -Wn 

dq 




From examination of equations (30), (32), (34-)> (35)> and (38), it 
becomes apparent that, for a given flight condition, calculation of the 
rate of change of A n with respect to five parameters (Bay, C nr , C n ^, 

C no , C 7 Q , and Cy q ') allows determination of the remaining partial 
P P P * 

derivatives, the amplitude coefficients associated with the lateral 
modes in the transient responses to unit yawing-moment, rolling-moment, 
or side-force inputs, and also the ratios which exist between rolling, 
yawing, and sideslipping in each lateral mode. 

Although the preceding analysis assumed unit moments or forces as 
inputs, the same type of approach can be used for arbitrary inputs with 
similar results. 
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The Slope Expression as a Possible Indication of Suitable 
Automatic Stabilizing Devices for Airplanes 

Flight tests and analysis of many of the present-day airplanes have 
i n dicated that the lateral or Dutch roll oscillation of these airplanes 
is very poorly damped. Usually, in order to improve this undesirable 
stability characteristic, artificial damping is supplied to the airplane 
by some type of automatic stabilization. However, the problem often 
exists that, although a specific stabilizer has proven to be satisfactory 
for one airplane, it has been found to be unsuitable for another. For 
instance, the damping characteristics of a certain airplane may be mark- 
edly improved by ah autopilot that supplies a yawing moment to the air- 
plane proportional to the rate of yaw but not much affected by an auto- 
pilot that supplies a yawing moment proportional to the rate of roll, 
and yet for another airplane the converse may be true. Also, this dif- 
ference In sensitivity of airplanes to various types of artificial 
dampers can exist even though they possess the same or nearly the same 
natural period and basic damping characteristics. Since the slopes or 
rates of change of the roots are a measure of the sensitivity of the 
damping of the airplane in the different modes of motion to various 
parameters, they should afford some insight into the conditions that 
might exist in order for some of the above situations to result. 

As an example, assume four airplanes which, for the same flight 
condition, have the same damping characteristics (same roots of charac- 
teristic equation) but differ in some of their masB and aerodynamic 
parameters. The characteristics of the assumed airplanes are given 
in the following table: 


Quantity 

Airplane 1 

Airplane 2 

Airplane 3 

Airplane 4 

b, ft 

35-3 

35-5 

35-5 

35-3 

V, ft/aee 

695-5 

695-5 

695-5 

695-5 

7 j de S 

0 

0 

0 

0 

c L 

0.24 

0.12 

0.12 

0.06 

t*b 

50 

50 

50 

50 

KX 2 

0.01485 

0.007425 

0.01485 

0.007425 

K Z 2 

0 . 0504 - 

0.1008 

0.0504 

0.1008 

K xz 

0 

0 

0 

0 

n, aeg 

0 

0 

0 

0 

Cjp, per radian . . 

-0.45 

-0.225 

-0.45 

-0.225 

C7 , per radian - 

0.04 

0.04 

0.08 

0.08 

t^, per radian 

-0.01 

-0.01 

-0.005 

-0.005 

Cn r , per radian 

-0.15 

-0.30 

-0.15 

-0.30 

Cyp, P er radian 

0 

0 

0 

0 

Cy r , per radian 

0 

0 

0 

0 

per radian ......... 

-0.58 

-0.58 

-0.58 

-0.58 

Cng; per radian 

0.12 

0.24 

0.12 

0.24 

Cjjj, per radian 

-0.11 

-0.11 

-0.22 

-0.22 
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Although the roots of the characteristic quartic for these airplanes are 
the same, their sensitivity to certain parameters is different. Examina- 
tion of the slopes should indicate then not only the parameters that are 
most effective in increasing the damping of the Dutch roll oscillation 
but also the reason that some parameters are_ more effective for one air- 
plane than they are for another. Since the quartics for the airplanes 
have the same roots, the coefficients of the quartics must also be 
identical. Therefore, if the slopes are to differ, the reason must be 
a difference in the partials of the coefficients of the stability quartic 
with respect to the various parameters. (See eq. (4).) 

Wow, assume that it is desired to equip each airplane with an 
auxiliary damping device in order to improve the Dutch roll damping. 

If the dampers are assumed to have no lags, the effect of the dampers 
can be considered merely As a modification of the stability parameters 
and the slopes will present a measure of the sensitivity of the Dutch 
roll damping to these parameters . 

The values for the rate of change of Dutch roll damping with respect 

to C n , C n , and C j for the four airplanes are given in the following 
“ P P 

table : 


Airplane 

da 

^Cn r 

da 

dCnp 

da 

aoi p 

1 

0.048 

-O.O 76 

O.OO 58 

2 

.024 

-.076 

.012 

3 

.048 

-.15 

.OO 58 

4 

.024 

-.15 

.012 


Assume, for example, that airplanes 1 and 3 equipped with a 
Cn r type yaw damper are satisfactorily stable. For airplanes 2 and 4, 

the same yaw danper would be only about 50 percent as effective as for 
airplanes 1 and 3* For all four airplanes it appears that, based on the 
slopes, a C n p damper would be more effective than a C nr damper, but, 

in selection of a suitable damper, consideration must be given to its 
effect on the other roots of the characteristic equation as well as to 
its effect on other factors such as the roll-to-yaw ratio. A Cn^ type 

damper, although generally effective in stabilizing the airplane Dutch 
roll, primarily redistributes the total damping of the system and, as a 
result, the damping- in- roll root A 2 is adversely affected. Hence, the 
value of which can be tolerated before this adverse effect becomes 

important is usually small and the increase in Dutch roll damping is 




NACA TN 313^ 


29 


somewhat limited. This adverse effect does not occur for the Cn r type 

damper and, consequently, a larger increase in Butch roll damping usually 
can he obtained by introducing large increments to Cn r • 


Nevertheless, the important point to be made here is that a dif- 
ference may exist in the sensitivity of their damping to different 
parameters even though the airplanes initially have identical damping 
characteristics. Futhermore, by use of the slope equation this differ- 
ence can be attributed to specific terms. For example, the reason air- 
planes 3 and 4 have a slope Sa/c£np double that of airplanes 1 and 2 


is that the values of the parameters C^ and Cj r for the former air- 
planes are exactly double these parameters for the latter airp lan es. 
Likewise the difference in the Cn r and C^p slopes can be attributed 

to specific terms. It is of interest to note that the difference in 
Sa/dCnr and Sa/SC n for the four airplanes will be predicted correctly 


by the simplified expressions (eqs. (24) and (26)). 


APPLICATION OF METHOD TO THREE AIRCRAFT AND 
DISCUSSION OF RESULTS 


Stability derivatives and mass and dimensional characteristics of 
three airplanes are given in table Il(a) . The nondimens ional roots of 
the quart ic equation (eq. (l)), the values of i/t-j^ for the aperiodic 
modes, and the values of 1 /t^^j and frequency of the oscillatory mode 
for each airplane are presented in table Il(b). 


The results obtained by evaluating the exact equation (eq. (4)) 
dA/dxq for the real roots and the complex roots for each airplane with 
respect to the ten parameters are given in table III. The rates of 
change of Dutch roll damping and frequency as calculated from equa- 
tions (17) to (22) are presented in table IV along with similar data 
from table 111 as a comparison of the results obtained by utilizing 
these simplified equations. Generally, the agreement Is very good. The 
slopes presented in tables IH and IV are nondimens ional slopes and If a 
comparison is made between the slopes of one airplane and those of 
another airplane, they must first be multiplied by the V/b ratio of 
the respective airplanes. These slopes are plotted in figures 1 to 4. 

In a lateral-stability analysis, T^yg :i - 8 commonly ns e d as 811 indi- 
cation of the degree of damping of an airplane. Here 1/T1/2 is chosen 
for plotting rather than Tp/2 since the former can be expressed as a 
linear function of Ax^ . 
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For the purpose of comparing the effect of increments in the param- 
eters on the lateral stability of the three airplanes and mainly as a 
comparison of the approximate changes , obtained by the methods of this 
paper, with the exact changes, figures 1 to 4 are presented. 

The straight lines in figures 1 to ^ indicate the approximate varia- 
tions of the ordinate due to various parameter changes, calculated by 
the methods of this paper. The symbols indicate check points which were 
taken at equal positive and negative increments for each parameter. 

These points show actual variations which were obtained by increasing 
or decreasing the specific parameter in the lateral equations of motion 
and solving the characteristic quartic (eq. (l)) for the new roots. The 
check points are presented in the following intervals, which are the 
same for each figure: between Axi /x± equal to ±100 percent for Cjp, 

Cl , and C nT J between Ax-t/xi equal to ±50 percent for C n „, Ci R , 
r r . p p Pq 

and Cy^J between Ax^/x^ equal to ±20 percent for Kx 0 and ; 

between AC^ equal to ± 0.05 per radian; and between At] equal to ±6°. 


Therefore, it is important to note that, with the exception of 
and T], the variations in l^T^/g 601 not plotted against 

equivalent incremental scales for the three airplanes since the basic 
parameter xi is generally different for each airplane. These particu- 
lar intervals were not chosen as limits or boundaries for the method but 
were used in an attempt to show the degree to. which the method could 
predict changes for. relatively large variations of many of the parameters. 


Figure 1 shows changes in the 



of the spiral root Aq of 


three airplanes due to increments in various airplane parameters. In 
figure l(a) at ACi^jCi^ equal to -100 percent, the actual variation 

has been plotted only for airplane B, since for airplanes A and C the 
real roots (Aq and A 2 ) have combined to form another oscillation. 


The spiral root of all the airplanes is most affected by those 
parameters which would be expected to cause the greatest changes as 
indicated from the approximate expression of the spiral root -E/D 
where E is the constant and D the coefficient of A in equation (l). 
One of the most important terms in the D coefficient is -(iCjpCnp. 

As seen from figure 1, and the parameters which constitute the 

E coefficient cause the most pronounced changes In the spiral root. 


The method appears to give good approximations to the changes In 
the l/T-jyg even for large variations of the parameters except for 
Cjp and C n p. For C the method predicted the changes fairly well 
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equal to ±4-0 percent. For Cn^ , it gave fair approxima- 


up to AC ip /C lp 

tlons up to ACn^y/Cnp equal to ±4-0 percent. 

Figure 2 shows changes in the 1/^/2 of the damping-in-roll 
root A 2 of the three airplanes due to increments in various airplane 
parameters. The method gave good agreement to the actual changes in 
1/^/2 for all variations of the parameters considered. In figure 2(a) 


at 


AC Zp /C 


equal to -100 percent, again the actual variation has been 


plotted only for airplane B. The reason for not plotting the actual 
variations for airplanes A and C was explained in the preceding discus- 
sion of figure 1. 


A good approximation to the damping-in-roll root is Ag 


Vb K X 2 

Therefore, from this expression, A 2 appears to be mainly sensitive to 
changes in Ci p and Ky^. Examination of figure 2 bears this out. 
Although Ky^. was not plotted, the variations of this parameter will 
almost be equal to the variations of Ky 0 ^. 

Figure 3 shows change In l^T^/g the Dutch roll oscillation of 

the three airplanes due to increments in various airplane parameters. 

The variation of Dutch roll damping with parameters other than Cj p 

is seen to be almost linear for each of the airplanes up to the largest 
changes considered for the parameters, and the agreement between the 
calculated damping and the damping predicted by the slope equations is 
seen to be very good. For airplane A, however, the effect on the spiral 
root Ap (see fig. l) and the Dutch roll damping (fig. 3(a)) of varying 
ACz p /C ip in the negative direction Is predicted by the slope equations 

only up to AC jpy/Cjp of approximately -4-0 percent. The calculated 

points indicate a reversal of the initial effect of Cj p on the Dutch 

roll damping beyond this value. The physical reason for this reversal 
Is not clear. Mathematically, however, the reversal can be shown to be 
the result of the relationship between the variation of the total damping 
of the system (b/a) and the variation of the sum of Ap and Ap with 
Cj . This line of reasoning leads to the conclusion that the damping pf 

£r 

the Dutch roll oscillation of =an airplane will be most sensitive to 
variations in Cj when 


d(-Ap - A 2 ) 

A 

d(B/A) 


T 
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or when these partials are equal or nearly equal and 


- ~Kq) 

r , 2 


is 


large. Although the reversal of the effect of Cjp on the Dutch roll 

damping of airplanes B and C is not apparent from figure 3(a); additional 
calculations indicated that this reversal does occur for larger varia- 
tions in than were included in the figure. 


Figure 4 shows changes in the frequency of the Dutch roll oscilla- 
tion of the three airplanes due to increments in various airplane 
parameters . The variations in frequency calculated from the slope 
equations are seen from these figures to agree almost identically with 
the actual variations for all of the parameters considered. A generally 
satisfactory approximation of the nondimens ional frequency of the Dutch 
roll oscillation is 


CD 


2Pb K Z 2 


K XZ C Ip N 

2p b Kx 2 Kz 2 


From the above expression the frequency appears to be mainly sensitive 
to Kz 2 ; Kx 2 , Kxz; and (and ^ which was not considered 

in this paper). Examination of figure 4 bears out the validity of the 
expression as an approximation to cd for the airplanes considered. The 
parameters Cn^ and Kzq 2 app.ear to have the greatest effect on cd 
but it should be pointed out that Kxz for the airplanes was generally 


small and therefore the term 


- K XZ 

^b K X £K Z 2 


contributed very little to cd 


for these airplanes. The general effect of this term can be seen from 
figure 4(d). As tj is assumed to be increased positively, the term 
becomes more positive and cd iB increased somewhat. 


CONCLUDING REMARKS 


Expressions have been presented in this paper which afford a means 
of calculating the rate of change of the roots of the lateral-stability 
quartic equation with respect to ten mass and aerodynamic parameters 
upon which these roots depend. Application of these expressions to 
three airplanes indicated that satisfactory results are obtained for 
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large variations in the ten parameters even though in the derivation of 
the expressions small variations in the parameters were assumed. An 
attempt was made to simplify certain of the expressions and the results 
indicated good agreement when compared with the results of the exact 
expressions . 

The expressions which define the rate of change of the character- 
istic stability roots with respect to the various parameters were shown 
to bear a definite relationship with the amplitude coefficients (and 
hence the ratios) of the lateral modes of motion subsequent to applied 
moments and forces. From these relationships calculation of the rate 
of change of 7^ with respect to five prescribed parameters allows 
determination of the remaining partial derivatives and these amplitude 
coefficients and ratios. 

The slope expressions should afford some insight into the types of 
automatic stabilization suitable for a given airplane since, if the 
dynamics of the automatic stabilizer are assumed to be negligible, the 
automatic stabilizer is effectively varying one or more of the mass or 
aerodynamic parameters of the airplane. 

The effect of a parameter change on a given airplane may be dif- 
ferent from the effect of the same parameter change on another airplane. 
The reasons for this difference can be seen from examination of the 
slope equation and, in some instances, can be attributed to particular 
parameters or combination of parameters. For example, the effect of 
the damping- in-yaw stability derivative Cn r appears to depend primarily 

on the moment of inertia in yaw, whereas the effect of the yawing-moment 
coefficient due to rolling velocity C n tends to be most effective for 

J. 

airplanes with low moment of inertia in roll or high effective dihedral. 
The relative effectiveness of C n ^ and C n ^ can be measured by the 

roll-to-yaw ratio ff \ |r of a given airplane, and for high values of 
0/i]r, Cnp is the most effective of the two. 

A method for calculating approximately the roots of the stability 
characteristic quartic is presented in the appendix and is shown to 
give very accurate results for the three airplanes considered. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 2, 1955- 
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APPENDIX A 

APPROXIMATE SOLUTION OF LATERAL STABILITY QUARTIC 


The period and damping of the various lateral modes of motion are 
determined from the linear and distinct roots of an equation of the 
form 


F(A) = AA4 + BA3 + CX2 + DX + E = 0 (Al) 

where the coefficients A, B, C, D, and E are functions of the sta- 
bility derivatives and mass parameters of the airplane. Rewriting, 
after division by A, gives 

F(A) - A 1 * - + B'A 5 + C'A 2 + D'A + E' =0 (A2) 


Assume that 


Hence , 


(A 2 + p^A + q^)(A^ + pgA + q.g) = F(A) 


B' 

c' 

D' 

e' 


Pi + P 2 

9-1 + 9-2 + ^ 1-^2 


PlO.2 + P29-1 
*1*2 


> 


(A3) 


Equation (A2) will almost alwayB have two real roots and a pair of con- 
jugate complex roots. One of the real roots Aq is associated with 
the inherently very poor spiral stability of airplanes and hence is 
very small. The remaining real root A 2 Is associated with the damping 
In roll of the airplane and Is of the order of the coefficient B ' . The 
real and imaginary parts of the complex roots determine the damping and 
period, respectively, of the normal lateral Dutch roll oscillation. If 
the assumption is made that the complex roots are obtained from the 
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quadratic A^ + p^A + q^_, the condition generally exists that qg « q^, 
and hence qg can be assumed to be zero in equations (A3)» Therefore, 


B' = p x + p 2 
c' “ q .1 + P 1 P 2 
D* = P 2 q. L 
E' ■ 0 


m 


From equations (Ph) the following expressions are obtained for p-^: 






(-A5) 


Pi 


_ ^(C* - q x ) 
D' 


Equating equations (A 5 ) yields 



Q-pCtlq - C ' ) 



Fl(q.l) = F 2 (qi) 


(A 6) 


Therefore, q^ can be determined by a simultaneous solution of equa- 
tions (A 5 ) . The functions F]_ and F 2 are illustrated by the following 
sketch: 
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Substitution of this value of q^_ into the first of equations (A5) 

gives 


Pi = 


C'fB'C' - D* ) 
(C'f + B’D’ 


(A9) 


or, since pi = -2a (where a is the real part of the complex root 
\ = a + ico), 


1 CCb'C* - D') 

2 (c 1 ) 2 + b’d' 


(A10) 


Also, the imaginary part of the complex root is related to q-^ By 


oi = (q x - a 2 ) 1 / 2 


(All) 


The two real roots may be obtained as follows : 
From equation (Ak) 


P2 


*1 


and from equation (A3) 


<*2 = 


ll 

*1 


(A12) 


Therefore, the quadratic giving the real roots may be written as 


7 ? + 2l A + ll 


^1 


*1 


0 
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Solving the above quadratic for A gives 


a - 4 { e ’ * ^’) 2 -^] 1/2 ) 


Then, the spiral root is given by the root with the plus sign 

r i /z\ 

\ ^ D 1 + J(D’) 2 - J (A13) 

and the damping-in-roll root is given by the root with the minus sign 

""N 

r 1/2 

A 2 = ^- <-D' - |(D') 2 - i^E*] r (Al4) 

Hence, the four roots of equations (l) may be found by equations (a 8), 
(AIO), (All), (A13), and (Al4). 

It should be reiterated that these results are contingent primarily 
upon the condition that « Q.q and if, upon determination of q_2 

(eq. (A12)), this condition is violated, the results are, of course, 
invalid, and resort must be made to other methods of determining the 
roots of the lateral-stability quartic,. 

The following table shows a comparison between the values of the 
nondimensional roots Obtained by using this approximate method and the 
exact roots for the three airplanes: 



Spiral root 

Damplng-in-roll root 

Dutch roll oscillation 

Airplane A 

Exact 

-O.OOO4107 

-0.13952 

-0.0094337 * 0.17271 

Approximate 

-.0004099 

-139*7 

-.009330 ± .17121 

Airplane B 

Exact 

-O.OOO 761 I 

-0.036142 

-0.00004245 ± 0.07091U 

Approximate 

-.0007646 

-03598 

-.OOOI 658 ± .071061 

Airplane C 

Exact 

- 0.00049 

- 0.15679 

- 0.00746 ± 0.156731 
-.00757 * .156861 

Approximate 

-.00049 

-.15656 
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TAB IE I— GENERAL FORMS OF THE PARTIAIS OF THE QUARTIC COEFFICIENTS 
FOR TEN DIFFERENT PARAMETERS - Continued 


% 

dVSGip 

0 

aB/ac lp 

-2ii b 2 Kz 2 

ac/ac lp 

^(^Z^p + I ^r) 

ao/ac Ip 

“ £ c n r CYp " HbCnp + ^ c n p CY r 

se / ac lp 

0 

% 

SA/ac Yp 

0 

^B/dCYp 

-W(kx%Z 2 - Kxz 2 ) 

ac/sc Yp 

R^X^r + " K XZ c *r " K XZ c n p ) 

ao/ac Yp 

^(Cn p Cz r - C^Cj^ 

SE/^Tp 

0 

Cnp 

dA/ac np 

0 

0 

Sc/SCnp 

(*b(jH*b% 2 + %Z C r p - Kx^Xr) 

3D/3Cnp 

-^b C i p + 2^-b c lfe + £ c i p c Y r ~ £ c Z r c Y p 

SE/SCnp 

- § C LCl r 
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TABI2 I— (3HERAL FORKS OT TSE PAKTIAIS C 9 TEE QUAETIC COSFICIEinS 
FOR TEH DIFFEttHT PARAMETERS - Concluded 


I ZkJhCl^ 

‘p 1 

c 1 

0 

as/acip 

0 

ac/acjg 

-ft^ft^z + Aa^p - K xz°rr) 

*>/acj p 

ftCbp - ^Cl-^Z 2 - ^ 2 npCy r + £ c n T 0'f p 

i &E/aci 0 

|°lPnr 

L : * 

aA/di^ 2 

Sft 5 ^ 2 

to/dl^o 2 

-2Mb S ^KZp 2 + Cnj. ccs^n + Ci p »ln 2 n + (cj r + C^sln n cob nj 

ac/ai^o 2 

ft [«« 2 n(Cn r 0T s + ^ftCnp - O^C np) + Bln 2 n(ci p C Y[) - Cy^) + 


am n «>» n(ci r OT ? + *ftCi p + c„ p cy b - Cn p c Xp - ci r oi p )| 

aD/iKXo 2 

-SubCL(ct„ am 2 n + Onj Btn 1 ) con t|) 


0 

V 2 

sa/^ 2 

<*bV 

Sb/^ 2 

-2ft 2 jay^Kj^ 2 *■ Cn r sin 2 !) + C ip coa s ri - (cj r + ft^aln i) cob 

ac/^Zo 2 

ft + ‘‘ftCj, s - <*,%) + coa 2 !)^!^ - C! s Ci b ) - 


Bln n 00 a 1 ( 01 ^^ + 4ftC ig + - C n pCy j - C Y r %3 

^D/dZz 0 2 

cos 2 tj - Cnp Bin t\ cob n) 

*/*% c 2 

0 

n 

dA/dvj 

0 

a8/a n 

£ft 2 {^xz(0 ip - C^) + (c, r + ^(k^ 2 - Ocob 2 i,] 

dC/&T| 

ft[2Klz(Cc-CT 3 + AftCnp - OjpCyp *0^ - C^) - (x^ 2 - + 


'‘ftClp + C^Oip - 0p s 0 Ip - Or r C lp )co. 2nj 

&)/&!* 

2ftC L [aCxzCjp + CjJsz/ - Kxo^oos 2^J 

2a/5q 

0 
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TAB IE II.- CHARACTERISTICS OF AIRPLANES CONSIDERED 

(a) Stability derivatives and mass and dimensional characteristics 

of three airplanes 







Airplane A 

Airplane B 

Airplane- C 

Wing: 

Span, ft . . 

• 




28.00 

25.00 

35-30 

Area, ft 2 . . 

. 




130.00 

175.00 

250.00 

Aspect ratio 

• 




6.00 

3.60 

^•97 

Sweepback, deg 

• 



5.05 

35-00 

0 

Dihedral, deg 


• 



0 

-3.00 

14-.00 

W j lb ■ • • • • 

W/S, Ib/ft® . . 

• 

• 



8,14.50.00 

9,214-5.00 

12,630.00 

• 




65.00 

52.80 

50.50 

Altitude, ft 

• 

• 

• 


30,000.00 

50,000.00 

30,000.00 

V, ft/sec . . . 

• 

• 

• 


797-00 

776.OO 

695.OO 

Cjj •••«•« 

• 

• 

• 


O.23 

O.I4-9 

0.21+ 

[1^ •■•••• 

• 

• 

• 


80.70 

182.00 

50.00 

T}, deg .... 

• 

• 

• 


-2.00 

0.82 

0 

KX 0 2 

• 

• 

• 


0.00962 

0.01557 

0.0114-85 

K Zq 2 ..... 

• 

• 

• 


0.05135 

O.I560 

0.05014- 

K X 2 

• 

• 

• 


O.OO967 

O.OI56 

O.OII4-85 

K% • • • • • 

■ 

• 

• 


0.0513 

O.156 

0.05014- 

%z 

» 


• 


-0.00145 

0.0020 

0 

CXp, per radian 

• 


• 


-0.14-0 

-0.33 

-O.I4.5 

C lr , per radian 

• 


• 


0.08 

0.23 

0.04 

Cnpj per radian 





-0.02 

-0.05 

-0.01 

C nr , per radian 


• 



-0.40 

-0.69 

-0.15 

Cy.p, per radian 


• 



0 

0 

0 

Cy r , per radian 


• 



0 

0 

0 

Cy 0 3 per radian 
P 

Cj^, per radian 


• 



-1.0 

-0.58 

-0.58 


• 



0.25 

0.25 

0.12 

Cjg, per radian 


• 



-0.126 

-0.18 

-0.11 

~b 9Ji *y • ■ • • • 


* 



0 

0 

0 



TABES II-- CHARACTERISTICS OF AIRPLANES CONSIDERED - Concluded 


£ 


(*) 


NondimenBional roots of the characteristic quartic equation, with the 
respective values of (<°' ) 0 for the three airplanes 



Honfl 1 mensional roots, A 

Damping factor, 

(^1/2),, 

Airplane A; 

Spiral root 

Damping- in- roll root . . . 
Dutch roll oscillation . . 

-0.0004107 
-0.13952 
-0.0094337 * 0.171271 

O.OI687 

5.7220 

0.3875 

Airplane B: 

Spiral root 

Damping-in- roll root . . . 
Dutch roll oscillation . . 

-O.OOO7611 
-0.036142 
-0.00004245 ± O.O709IH 

0.03409 

1.6190 

0.001901 

Airplane C: 

Spiral root 

Damping- in- roll root . . . 
IXitch roll oscillation . . 

-0.00049 
-0.15679 
-0.00746 ± 0.156731 

0.01393 

4.4580 

0.2121 


Frequency, 

(®’) 0 


4.875 


2.201 


5.088 


1 ■> 
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TABIE HI.- HOKDIMEHSIOHAL SLOPES FOR THREE -AIRPLAKES WITH 
RESPECT TO TEH PARAMETERS 



Slopes 

x i 





Airplane A 

% 

-O.OOIO 

0.32 

0.0024- 

-O.OO9O 

Cir 

.0034 

.0026 

.0015 

.0040 

Cn r 

.0017 

.00088 

.029 

.0023 

Cn p 

-.0004-9 

.11 

-.048 

-.041 

% 

.0026 

.030 

-.016 

■ 54 

Clp 

.0052 

.088 

-.047 

.020 

% 

n 

-.0000019 

.00018 

.0030 

-.00013 

.000017 

.16 

-.088 

.098 

W 

.00021 

-.11 

.29 

-1.68 

KXo 2 

-.00012 

1A.47 

-57 

-.076 


Airplane B 



-0.0017 

0.077 

0.0062 

-0.0015 

Ci r 

.0029 

-.0023 

-.00086 

.00057 

c n r 

.0021 

-.0013 

.0040 

.00072 

% 

-.0012 

.044 

-.022 

1 

» 

O 

G 

% 

°lp 

.0048 

.018 

-.011 

.12 

.OO67 

.031 

-.015 

O 

• 

1 

% 

-.0000019 

.000094 

.0013 

.0000060 

n 

.00013 

.16 

-.089 

.16 


.0012 

-.019 

.028 

-.19 

K *o 

-.00094 

2.022 

-.068 

-.34 

Airplane C 

Ci P 

- 0.0010 

0.33 

O.OO58 

-0.00060 

Cl r 

.0050 

-.0016 

-.0017 

.0019 

Cn r 

.0046 

-.00074 

.048 

.0028 

% 

-.00094 

•15 

-.076 

-.077 

Cnp 

.0055 

.021 

-.013 

.62 

C h 

.0060 

.046 

-.026 

-.022 

-.0000028 

.00016 

.0049 

-.00013 

n 

.000014 

•17 

-.088 

cO 

0 

« 

V 

.00045 

-.023 

.16 

-1.49 

K Xo 2 

-.00010 

10.22 

-.010 

-.18 
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TABLE IV.- RATE OP CHARGE OF DUTCH ROLL DAMPING AND FREQUENCY OBTAINED 
FROM SIMPLIFIED EXPRESSIONS COMPARED WITH THE EXACT 


x i 

da/ 5 xi 

Stn/dx 1 

Simplified 

i 

Exact 

Simplified 

Exact 



Airplane A 


% 

0.0031 

0.0024 

-0.0086 

-0.0090 

ci r 

.0016 

.0015 

.0037 

.0040 

°nr 

.029 

.029 

.0017 

.0023 

c n p 

-.045 

-.048 

-.042 

-.041 

s 

-.014 

-.016 

• 3 4 

• 34 

°ip 

-.045 

-.047 

.019 

.020 

% 

.0030 

.0030 

-.00016 

-.00013 

n 

-.074 

-.088 

.098 

.098 

KZ 0 2 

• 29 

• 29 

-I.69 

-1.68 

%o 2 

-.54 

-•57 

-.10 

-.076 



Airplane B 




Cl P 

0.0062 

0.0062 

-0.0016 

-O.OO15 

C l r 

-.OOO85 

-.00086 

.00058 

.OOO57 

Cn r 

.0040 

.0040 

.00067 

.00072 

Cnp 

-.022 

-.022 

-.012 

-.013 

°np 

-.0099 

-.011 

.12 

.12 

cz p 

-.019 

-.019 

-.027 

-.024 

°Yp 

.0013 

.0013 

-.000009 

.000006 

n 

-.O87 

-.089 

.16 

.16 

% 2 

.026 

.028 

-.19 

-19 

*0® 

-.071 

-.068 

-•33 

-.34 

Airplane C 

°*P 

O.OO58 

0.0058 

-0.00021 

-0.00060 

Cl r 

-.0018 

-.0017 

.0018 . 

.0019 

C n r 

.048 

.048 

.0024 

.0028 

% 

-.073 

-.076 

-.077 

-.077 

c np 

-.011 

-.013 

.62 

.62 

°ip 

-.025 

-.026 

-.022 ... 

-.022 


.0049 

.0049 

-.00015 

-.00013 


-.083 

-.088 

.083 

.083 

V 

• 17 

.16 

- 1.49 

- 1.49 

K Xo 2 

-.013 

-.010 

-.17 

-.18 
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(a) 



and Cj r » 


Figure i.- Changes in the 



of the spiral root of three airplanes 


due to increments in various airplane parameters . 
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Airplane Approximate variation Actual variation 
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(b) and 

Figure 1.- Continued, 
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Airplane Approximate variation Actual variation 



(e) K^ 2 and K^ 2 . 


Figure 1.- Concluded 
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Airp lane Approximate variation Actual variation 

A - r\ v n 
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(a) C? and C 7 . 
ip l r 


Figure 2 .- Changes in the 



of the damping-in-roll root of three 


airplanes due to increments in various airplane parameters . 
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Figxire 2.- Continued. 
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* leare a- Continued 
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Figure 3.- 

Continued. 
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(d) Oy and tj. 
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figure 5.- Continued 





Airplane Approximate variation Actual variation 



airplanes due to increments in various airplane parameters. 




Figure 4 -- Continued. 
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Figure k.- Continued. 
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Figure lj-.- Continued. 


1t£t£ Ni VOVN 





